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We propose and develop to some extent a novel approach, which allows us to effectively 
describe, for relativistic heavy-ion collisions, the empirically observed deviation from unity 
of the intercept A (i.e. the measured value corresponding to zero relative momentum p of 
two registered identical pions or kaons) of the two-particle correlation function C(p,K). 
The approach is based on the use of two versions of the so-called (/-deformed oscillators 
and the corresponding picture of ideal gases of g-bosons. By these techniques the inter- 
cept A is put into direct correspondence with the deformation parameter q. For fixed 
deformation strength, the model predicts dependence of the intercept A on the pion pair 
mean momentum K. 



X 



1 Introduction 

The hadron matter under intense conditions of high temperatures and densities has been ex- 
tensively studied with the use of relativistic heavy- ion collisions (RHIC). The insight into the 
extreme matter with experimentally controlled initial energies serves as an examination of ex- 
isting models and theories. On the other hand, RHIC promises to be a laboratory where a 
search for new physics that extends beyond current imaginations can be made. 

The models and approaches that are used to describe the processes occurring in the reaction 
region are examined by comparing provided predictions with experimental data on single-, two- 
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and many-particle momentum spectra, which contain information on the source at the early 
stage (photons, dileptons) and at the stage of so-called " freeze- ouf (hadron spectra). Two- 
particle correlations encapsulate information about the space-time structure and dynamics of 
the emitting source [[]]]- 0. Usually, consideration of the correlations that occur in relativistic 
heavy ion-collisions assumes that: (i) the particles are emitted independently (or the source is 
completely chaotic), and (ii) finite multiplicity corrections can be neglected. Then, correlations 
reflect a) the effects from symmetrization (antisymmetrization) of the amplitude to detect 
identical particles with certain momenta, and b) the effects that are generated by the final state 
interactions of the detected particles between them and with the source. At first sight, the final 
state interactions (FSI) can be regarded as a contamination of "pure" particle correlations. 
However, it should be noted that the FSI depend on the structure of the emitting source and 
thus provide information about source dynamics as well [|J . Discussion of the latter is beyond 
the scope of the present paper. 

The nominal quantity expressing the correlation function in terms of experimental distri- 
butions is 

c ^™ = !$£&' (1) 

where Pi (k) = Ed 3 N/d 3 k and P 2 (k a , k&) = E a E b d 6 N/(d 3 k a d 3 k b ) are single- and two-particle 
cross-sections, k a and k& being on-shell asymptotic momenta. 

In the absence of FSI, for a chaotic source, the correlation function can be expressed as (see 
Appendix): 

\fd 4 Xe l P- x S(X,K)\ 2 

C(p,K) = 1 + cosa f \ I V > ( 2 ) 

fd 4 XS(x,K + 2) f<PYS[Y,K-l) 

where 4-momenta K and p defined as 

K = -(k a + k b ), p = k a -k b . (3) 

The source function S(x, K) (single-particle Wigner density) is defined by emitted single- 
particle states ^ 7 (£,x) at freeze-out times: 

S(Y, K) = J d 4 y e iK -y £ p lY ^ (Y + f ) ^ {Y - f) , (4) 

7:7' 

where the summation (averaging) is taken over the set of all quantum numbers {7} carried 
by the particle just before it is frozen out. The source freeze-out density matrix p 7 y is the 
weight factor of such an averaging and depends on the particular model of source, for instance, 
thermal density operator is widely exploited. 

For the system of identical particles that we are going to consider, the two-particle wave 
function appears to be a symmetrized (antisymmetrized) construction of single-particle states 
(chaoticity assumption), which reads 

^ 7 a76 ( X "> x b, Q = ~j= [VV (*a, *a) ^7b ( X &> *a) + C^a ( x b> *o) ^76 ( X a, <a)] , (5) 



where a = for identical bosons, and a = ir for identical fermions. Note that the function (|5]) 
is taken at freeze-out times (about translation from detector times t — > oo to emission times, 
see Appendix). 

From now on we shall refer our consideration of two-particle correlations to identical bosons 
(pions, kaons, etc.). As follows from Eq. (|2|), the boson correlation function should approach the 
exact value two as the relative momentum approaches zero. But as it was observed, from the 
very first experimental data and up to the most recent experiments, the measured correlation 
function never reaches this value at p = 0. To remove this discrepancy, the correlation function 
of identical bosons is always taken in the form 

C(p,K) = l + Xf(p,K), (6) 

where A is drawn from an experimental fit to the data, usually in the range A = 0.4 - 0.9; 
f(p,K) is commonly taken as a Gaussian function (in any case, /(p = 0, K) = 1). The 
deviation of A from unity in RHIC can be explained by the production of secondary pions from 
resonance decays which are outside the fireball. The presence of long-lived resonances results 
in an increase of the measured source size and life-times J5], [| . 

We are now coming to the key idea of our paper. Confining ourselves to formula (Q) to 
try to explain experimental data, it is then straightforward to put into correspondence the 
parameter A with the angle a, so as to get, by means of cos a, the right reduction factor A. 
Actually, in Eq. (5) one can take the phase factor e~ ta in place of the factor e ta . However, 
by simple algebra the two-particle amplitude (52) can be reduced to the form which results 
in the same two-particle probability as the former one. Indeed, the correlation function (0), 
which is a measurable quantity, is obviously symmetric with respect to a — > —a. It turns out 
that an effective symmetrization of a two-particle wave function in heavy-ion collisions exhibits 
similar features to what one encounters in the description of the Aharonov-Bohm effect, or in 
the physics of anyons. This means that the two-particle wave function of a boson pair released 
from a dense and hot environment effectively acquires an additional phase. Hence, the drawn 
phenomenon can be ascribed to the properties of the medium formed in RHIC, which, as we 
see, exhibits some non-standard QFT behaviour through the considered correlation functions. 
So, adopting as a driving idea the fact that the correlation function approaches 1 + A when the 
two-boson relative momentum approaches zero, we will attempt to construct an effective model 
capable to mimic the real physical picture. To perform this, we shall use as our basic tool 
the so-called (algebra of) q-deformed commutation relations, or techniques of q-boson statistics, 
which certainly can be put in connection with the symmetrization rules. 

The deformation parameter q is viewed as an effective (not universal) parameter which effi- 
ciently encapsulates most essential features of complicated dynamics of the system under study. 
In many cases, usage of appropriate g-algebra allows one to reduce the treatment of complex 
system of interacting particles to consideration of a system of non-interacting ones at the price 
of complication (deformation) of the commutation relations. As an example let us mention the 
application of g-deformed algebras to description of rotational spectra of superdeformed nuclei 
- here q has different value for each nucleus. 

It is worth noting that in the context of hadron theory, g-deformed algebras (or quantum 
algebras) were also already applied. Such a usage proved to yield a significantly improved 



description of hadron characteristics, both regarding hadron scattering ||-[1C] - nonlinearity 



of Regge trajectories - and in the sector of such static properties as hadron masses and mass 



sum rules [[Ll], [12| . 



In what follows, we shall exploit, for the system of pions or kaons, the (ideal) g-Bose gas 
picture based on two concrete versions of g-bosons. From the viewpoint of direct physical mean- 
ing and/or explanation of the true origin of the g-deformation in the considered phenomenon, 
these versions differ from each other, first of all in the question of whether q must be real or 
can also take such complex values as a pure phase. Note that, at this stage, we do not go into 
details concerning the diversity of other, than the above-mentioned, "microscopical" reasons 
(certainly, not completely unrelated) for the appearance of g-deformed statistics. Suffice it to 
mention that the composite nature of the particles (pseudoscalar mesons) under study may as 



well result [13, 1 [\ in the g-deformed structures linked to the real deformation parameter g. 



2 The two versions of q-bosons 

In this section we give a brief sketch of main features of the two versions (type "A" and type 
"B") of multimode g-oscillators, which will be used in subsequent treatment. 

Type A 

The g-oscillators of this type are defined by means of the relations [[HJ 

[a h %] = [a], a]] = , [A/*, a,-] = S^ttj , [A/i, a]] = S^a] , [A/i,A/}] = , 

aia] -q Sij a]ai = 8% . (7) 

Note that, if z ^ j, this system of independent g-oscillators differs essentially from quons [HJ 
whose different modes are non-commuting (i.e. g-commuting). 

From the vacuum state given by Oj|0, 0, . . .) = for all i, the state vectors 

Kna Bi,-)s /r ... * ===(o}r(c4r---(ot)"«-..|0 > 0,...) (8) 

are constructed as usual, so that 

a\\n 1 ,...,n i ,...) = ^/\n~+i}\n 1 ,. . . ,n { + 1, . . .) , (9) 



Oi|ni, ...,7ii, ...) = V[ni]|ni,...,nj - 1,...) , (10) 

A/i|ni,...,7ij, ...) = m\ni,...,rii, . . .) . (11) 



Here the bracketed notation 
r i * - <f 



along with [r]! = [1][2] ■ ■ • [r - l][r] , [0]! = 1, (12) 



is used. The g-bracket [A] for an operator A is understood as a formal expression (formal 
series). At q — > 1, from [r] and [A] we recover r and A, thus returning to the formulas for the 
standard bosonic oscillator. In what follows it will be assumed that 

-i <<?<!; (13) 

for each such value of the deformation parameter q, the operators a] , Oj are mutual conjugates. 

In the generic case where q ^ 1, the bilinear ajaj does not equal the number operator Mi 
(as this is true for usual bosonic oscillators) but, instead, 

a\ ai = [Mi\. (14) 

The inverse of the latter relation is given by the formula [ |TTjj : 



^ = Et^(«!)x, (15) 

s=l 

expressing the number operator as a (formal) series of creation and annihilation operators. 

TypeB 



The g-oscillators of the second type are defined through the relations |L8|, [19| : 

[k, b 3 ] = [6j, b)\ = , [N t , b,j\ = -5 tJ b j , [Ni, b]] = Sijb] , [Ni, NA = , 

kb) - q 5 *b)k = 6 ijq - N ' , hb) - q-^btbi = 5 ijq N i . (16) 

Again we have 

h% = [Ni] (17) 

where the notation for the (/-bracket this time means: 

V — V 

q — q 



q-q 1 



(18) 



Formulas completely analogous to Eqs. (8)— (11) are valid also for the operators bi, bj if, instead 

of (12), we now use the definition (18) for g-brackets. Clearly, the equality b\bi = N holds only 
in the "no- deformation" limit of q = 1. For consistency of the conjugation, it is required that 
either q is real or 

q = exp(i6) , < 6 < 2vr . (19) 

In what follows we will consider the type "A" as well as the type "B" oscillators. For the 
type "B" oscillators the exponent form (|T9|) will be adopted (compare with the symmetrization 
phase a in Eq. (BI)). 



3 Statistical (/-distributions 

For the dynamical multi-pion or multi-kaon system, we consider the model of ideal gas of 



g-bosons (IQBG) taking the free (non-interacting) Hamiltonian in the form |2(| |21 



H = J2^iM, (20) 



where oot = ^Jra 2 + k^ , A/i is defined as above, and subscript i labels energy eigenvalues. It 
should be emphasized that among a large variety of possible choices of Hamiltonians, this is 
the unique truly non-interacting one, which possesses an additive spectrum. From now on, we 
assume that 3-momenta of particles take their values from a discrete set (i.e. the system is 
contained in a large finite box of volume ~ L 3 ). 

As usual, basic statistical properties are obtained by evaluating thermal averages such as 

where j3 = 1/T and the Boltzmann constant is set equal to 1. The averaging here is taken with 
respect to the chosen Hamiltonian (20). 

It is an easy task to calculate the quantity (q ), and to obtain 

(q*) = % . (21) 

w ' e#* -q 

From this we find the distribution function (recall that q is from the interval — 1 < q < 1): 



<«!«*> = -J— ■ (22) 



In the no-deformation limit q — > 1, this reduces to the Planck-Bose-Einstein distribution, as it 
should, since at q = 1 we return to the standard system of bosonic commutation relations. 

At q = — 1 and q = 0, the distributions we get coincide respectively with Fermi-Dirac and 
Maxwell-Boltzmann ones. It should be emphasized that this coincidence is rather formal: the 
defining relations (7) at q = — 1 or q = differ from those for the system of fermions or the non- 
quantal (classical) system. The formal coincidence of Eq. (22) at q = — 1 with a Fermi-Dirac 
distribution can be interpreted |22| as due to the impenetrability (the hard-core property) of 



such bosons. The difference with the system of genuine fermions lies in commuting (versus 
fermionic anticommuting) of non- coinciding modes at q = —1; see (7). 

Let us now turn to the type "B" g-bosons. The Hamiltonian is chosen again as that of 
IQBG with the number operator defined in (16) and (17), i.e., 



H = ^u i N l . (23) 



Calculation of (q ±Ni ) yields 



e 



fat 



(^) = -T- -TT • (24) 



e /3u>i _ q±i 



From the relation 



(b]bi) 



-Ni 



o/3u>i 



or, equivalently, (b]bi) = (alai)BE{q Ni )(q Ni ), we obtain the formula for the g-deformed distri- 
bution function (taking into account that q + g _1 = [2] = 2cos#): 



(b\h) 



oP^i 



o2f3u>i 



2cos(0)e#* + 1 ' 



(25) 



Note that, although the deformation parameter q is chosen in a particular complex form, see 
(19), the explicit expression for the g-distribution function turns out to be real, owing to its 
specific dependence on q through the combination q + g _1 . 

The shape of the function f(k) = (tfb)(k) from (25) corresponding to the gas of pions 
modelled by IQBG is pictureQin Fig. 1 (curve II). For comparison, the standard Bose-Einstein 
distribution function (curve I) and the classical Maxwell-Boltzmann one (curve III) are also 
presented in the same figure. As is clearly seen, the g-deformed distribution function lies 
completely in between the other two curves, thus demonstrating^ that the deviation of the 
g-distribution (25) from the quantum Bose-Einstein distribution goes in the "right direction", 
towards the classical Maxwell-Boltzmann one. 
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Figure 1: The g-distribution function (25) versus momentum (curve II), in comparison with 
the quantum Bose-Einstein (curve I) and classical Maxwell-Boltzmann (curve III) distributions. 
The inputs are: T = 120 MeV, m = m^ ; curve II corresponds to the deformation angle 9 = 24°. 



4 The (isotriplet-averaged) pion mass m(7r ±,c) ) = 139.57 MeV and the temperature T = 120 MeV are taken 
as inputs. The deformation value, encoded in cos#, is fixed to be 8 = 24°. 

5 The same is true also for the (apparently more simple) g-distribution (22) of the type "A" g-bosons. 



Analogous curves for g-distribution functions, with similar properties, can be given for other 
fixed data. Let us remark that, in the case of kaons, because of their larger mass and higher 
empirical value of the intercept A ~ 0.88 (which corresponds to the smaller deformation in our 
model), such a curve should lie significantly closer to that of the Bose-Einstein distribution. 

It is worth noting that the g-distribution functions (22) and (25) already appeared (in the 



context of thermal field theory) in |2l| . We give them here for completeness of exposition, since 



they are connected with new results to be described in the next section. 

4 Two-particle correlations of ^-bosons 

Let us now turn to the issue of two-particle correlations. From the easily verifiable identity 

a\a]a k ai - q~ hk ~ kl a}-a k a { a\ = [a\, a]]a k ai + aj[aj, a k ] q ~s lk ai + q~^ k a]a k [a\, a{\ q -s u , 
by taking thermal averages, we find 

{a\a)a kai ) = ^3— J -^-—((a]a l ) (a\a k ) + q~ 5 ^ (a]a k ) (a]^)) . 

With coinciding modes^] , this leads to the formula 

t I 1 + q 

(a\a\aiai) = — ~ -— — . (26) 

N l % ' (e#* - q)(e' 3ui - q 2 ) 

From the last relation and distribution (22), the ratio under question follows: 

~ _ (tJoJojOj) _ (1 + q)(e^ - q) 
l ~ (al^) 2 ~ e^-q* ' ^ 

This constitutes one of our main results. For convenience, let us set 

A = 1 + A with A = q— a . (28) 

e^ - q 2 

The quantity A can be directly confronted with empirical data. Note that in the non-deformed 
limit q —* 1 the value Abe = 1 5 proper for Bose-Einstein statistics, is correctly reproduced from 
Eq. (28). This, obviously, corresponds to the Bose-Einstein distribution contained in (22) if 

g-> 1. 

At q = —1, Eqs. (22) and (28) formally coincide with the Fermi-Dirac distribution and 

the value Afd — — 1 proper for the Fermi-Dirac statistics respectively, although the defining 
relations are not exactly those of the fermionic system (rather, the hard core bosons |22]j ). 
Finally, at q = we get A = 0, which coincides with the analogous fact for the case of purely 
classical description (complete absence of quantum effects due to identical particles). The three 
different cases are clearly seen in Fig. 2 as the only three points where all the different curves 



Recall that K = k(k a + kt), p = k — kt; in the case p = 0, K = k a = k& 




Figure 2: Intercept A versus deformation parameter q, as given by Eq. (28). The curves A, 
B, C and D correspond to the values wa = 0.1, wb = 0.35, wc = 0.8 and Wd = 5.0 of the 
dimensionless variable w = /3u. 



(the continuum parametrized by to = /3u) merge and, thus, the dependence on momentum 
and/or temperature disappears. From the continuum of curves, there exists a unique limiting 
(asymptotic) one A = 1 + q (or A = q), which corresponds to the limit w — ► oo (i.e. to zero 
temperature or infinite momentum). Conversely, for very large temperature such that w — » 0, 
the curve goes over into the step-shaped function 



A 



with the constant A = for each fixed value of q except for the endpoints q — 1 and q — — 1. 

We find now the formula describing two-particle correlations (at identical momenta), which 
corresponds to Biedenharn-Macfarlane g-oscillators, see Eq. (16). From the relation 

(b\b\bA) ~ q 2 (b\bAbl) = -{b\b iq Ni ){l + q 2 ) , 
which is valid just for coinciding modes (i.e. equal momenta), we immediately get 




(b\b\bA) 



l + g 2 



«2g/3ui 



-<&w*> 



The thermal average in the r.h.s. can be easily evaluated to yield (b\biq Ni ) = q/(e /3u ' i — q 2 ) 



Taking this into account, we find the expression for two-particle distribution, namely 

(llw = e ^_ 2 ^)^- + i ■ (29) 

From that, the desired formula for (the intercept of) two-particle correlations finally results in: 

I = x , , (&MW 2cosfl(t t + l-cosfl) 2 
l ~ t+ «&JW~ «? + 2(1 " cos**)* ' (30) 

where tj = cosh(/9cjj) — 1. Note that both (29) and (30) are real functions (as they should) 
since, like (25), they both depend on the complex g-parameter (19) through the combination 
\{q + q~ v ) = cos9. 

In the rest of this section, we extract some useful information contained in Eq. (30). Solving 
this equation with respect to cos 9 at a fixed value A = A, we obtain the deformatin angle as 
the function: 9 = 9(X, K, T, m). 

Figure 3 is used to illustrate main properties of the intercept (correlation strength) A treated 
from the standpoint of g-deformation, that is, on the base of Eq. (30). First, let us note that the 
continuum of curves A = A(cos#) parametrized by w = f3uo divides into three different classes 
( "subcontinua" ) given by three intervals of the parameter: (i) < w < Wq, (ii) wo < w < w' , 
and (iii) w' Q < w < oo. Here, the two "critical" values w = wb — 0.481 and w' Q = w^ — 0.696 
are singled out (curves B and D respectively). The curves A, C, and E are typical representatives 
of the classes (i),(ii) and (iii). All the curves from classes (i), (ii) possess two extrema, the 
minimum being to the left of the maximum (the curve D is unique since its extrema degenerate, 
coinciding with the point of inflection). This fact enables us to define naturally "the range of 
small deformations" - the interval / sma ii for the variable 9: from 9 = (no deformation) to the 
value yielding minimal A, A m i n ~ 0.33, implied by the "critical" value wo = wd- It is seen that 
on the interval i sma ii the intercept A monotonically decreases with increasing of 9 (or 1 — cos 9), 
the strength of deformation. On the contrary, each curve from the third class is monotonic 
as a whole and, thus, there is no criterion (no peculiar point), which would naturally separate 
"small" deformations from "large" ones. 

In this paper we deal only with classes (ii) and (iii) since all their curves at 9 ^ lie below 
the straight line A = 2 - the largest possible correlation attainable in the Bose-Einstein case 
(note that curve B contains, besides 9 = 0, just a single point at a certain value of 1 — cos#, 
where the value A = 2 is also attained). Moreover, for the (ii)-type curves, we restrict ourselves 
to /small, ignoring "moderately large" deformations (between min. and max.), for which the 
behaviour of A = X(9) is opposite to that for I Bma a, as well as very large ones (to the right of 
the maximum). Unlike these two regular classes, the class (i) consists of "irregular" curves: 
for each such curve there exist g-deformations that generate correlation strengths exceeding the 
maximal possible one A = 2 . Therefore, we discard the class (i), at least at this stage. 

Finally, let us discuss special values of the physical variables T, |K|, which provide the 
peculiar values wq = wb — 0.481 and w' Q = wd — 0.696 (recall that w = y/m 2 + K 2 /T). With 
m(7r ± '°) = 139.57 MeV and lowest mean momentum of the pion pair fixed to be |K| = 0, we 
get two bounded from below values for the temperature: T = 290.0 MeV and Tq = 200.5 MeV. 

10 




Figure 3: Intercept A versus deformation given by cos8, see Eq. (30). The curves A,B,C,D 
and E correspond to the values wa = 0.3, wb = 0.481, wc = 0.58, wd = 0.696 and we = 2.0 
of the dimensionless variable w = f3u. 

It is interesting to compare these data with that for the typical curve from class (iii). Namely, 
at w — We = 2.0 (the curve E) for pions of this same lowest momentum we get the limiting 
temperature T E = 69.8 MeV. 



Discussion and outlook 



The main purpose of theoretical approaches to RHIC is to find an adequate description for the 
non-equilibrium state formed during the collision. On this way, the g-boson techniques enables 
us to treat the non- stationary hot and dense matter effectively as a "noninteracting ideal gas" . 
To deal with g-bosons, it is necessary to determine the g-parameter that corresponds to the 
actual state of the hot medium. We propose a way of extracting from the two-particle correla- 
tions a useful information concerning q, and develop an effective picture of the two-pion (-kaon) 
spectra in RHIC. According to our results, the measured deviation from unity of the intercept 
A is interpreted as the manifestation of g-boson properties of the pion system created in RHIC. 

From Eqs. (28) and (30), one can express the quantities encoding the deformation as: 
q = q(X, K, T, to) and 6 = 6(X, K, T, to). In the characteristic limits of temperature 
and momentum, our model exhibits a remarkable feature, valid for both types of exploited 
g-oscillators: 

• For very low temperature at fixed momenta, or very large momenta at fixed temperature 
(i.e. at w — > oo, compare with curve D in Fig. 2 and curve E in Fig. 3), we come to the equality 



q = \ 



(T^0 



or 



IKI 



oo 



(31) 
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with the type "A " oscillators, and to the equality 

2cos# = A + l 



(T^O 



or 



IKI 



oo) 



(32) 



for the type "B" oscillators. This implies the unified direct connection A <-» g, namely A = [2], 
for both types "A" and "B" . 

On the other hand, finite temperature and momenta become non-trivially involved (espe- 
cially in the case of type "B" g-bosons) in the relation between A and the deformation parameter 
q, see Eqs. (28), (30) as well as Figs. 2, 3. 

• Equation (30) implies dependence of the intercept A on the pair mean momentum K for 
fixed values of the deformation angle 9, temperature and particle mass, since A = X(9, K, T, m). 
In Fig. 4, we present this dependence exemplifying it, for T = 120 MeV and m = m n , with 
four curves which correspond to the values 4°, 9°, 15° and 24° of the deformation angle 9. Each 
curve tends to its own asymptote given by Eq. (32). The apparent variability of A with varying 
|K| (for a fixed deformation) is one of the consequences of our model, and can be viewed both 
as an interesting prediction and as a check-point for this approach. 




Figure 4: Dependence of the intercept A on the pion momentum |K| for some values of 9. The 
temperature is fixed to be T = 120 MeV; m = m n . 



In p 



it was argued that the quon-based description applied to pions participating in the 
decay K^ — > 7r + 7r - , can give a complete account of this process (instead of CP- violation, as 
commonly accepted cause). Moreover, this sets a bound on the "strength" of the deformation, 
namely 1 — q < 10~ 6 . On the contrary, within the proposed approach to multiparticle correla- 
tions (as occurred in the HBT interferometry) based on viewing the identical pions as g-bosons, 
the two-pion correlations imply the employment of significantly more developed (strength of) 
g-deformation. Namely, in our implementation of the g-Bose gas picture, the empirical data 
are interpreted as corresponding to the strength ~ 0.2 0.6 of g-deformation. 
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The approach proposed in this paper, we hope, opens interesting new perspectives for further 
research in this direction. Three-particle correlations as well as fermion-fermion(-fermion) ones 
are the topics to be studied next using the developed "g-techniques" . 
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APPENDIX 

In the appendix we consider the two-particle quantum statistical correlations when the final 
state interactions of the detected particles are neglected. This phenomenon is visualized most 
transparently on the basis of standard quantum mechanics in the non-relativistic case. Mean- 
while, using the presented scheme, the relativistic picture can also be considered. Moreover, 
this approach allows us to take into consideration the final state interactions as well ||]. 

The probability to register two particles, with definite asymptotic momenta p a and p b , which 
are created in the relativistic heavy-ion collisions is usually compared with the probability to 
register two particles of the same momentum independently. For that reason, we begin by 
considering the single-particle spectrum. 

1. Single-particle cross section 

Let us consider the single-particle state ^7 of a particle emitted by the source. Its propagation 
to the detector is governed by the Schrodinger equation 

^ 7 (x,t) = ~ ^^ ~ = l_ v2 ^ 

at 2m 

The index 7 denotes a complete set of 1-particle quantum numbers. (In a basis of wave packets, 
these could contain the centers X of the wave packets of the particles at their freeze-out times 
t.) Equation ([33]) is solved by 

^ 7 (x,t,t ) = e-^( 4 - 4o ^ 7 (x,t ), (34) 

in terms of the single-particle wave function at some initial time to- We will assume that the 
detector measures asymptotic momentum eigenstates, i.e. that it acts by projecting the emitted 
1-particle state onto 

(/.° ut (x,t) = e i(p - x ^ (p)<) , (35) 
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where w(p) = p 2 /2m is the energy of the particle. The measured single-particle momentum 
amplitude is then 

Ay(p, t ) = lim / d 3 x 0° ut '*(x, t) V 7 (x, t, t ) • (36) 

Using the time evolution equation fl34"|), this can be expressed in terms of the emitted single- 
particle wave function ip 1 at earlier times as 



AJp,t ) = lim [ d 3 x \e- ih ^ t0 -V ^{x, t) 

t^ooj L 



V> 7 (x,£ ) 



= /"rf 3 xe^ (p)f °- 4p - x ^ 7 (x,t ). (37) 

This expression means that the momentum amplitude, as it should be, is an on-shell Fourier 
transformation of the emitted wave at emission times. It is worth noting that this is not the 
case when, after emission, a particle is subject to final state interactions. 

The single-particle probability is obtained by averaging ($% and its complex conjugate using 
the density matrix defining the source. This density matrix is characterized by a probability 
distribution for the single-particle quantum numbers 7. We write 



Pi(p) = ^|A 7 (p,t )| 2 ) = J2 dt dt f p^,A y ,(p,t' ) A*(p,t ) 

77' 



(3t 



where the anglular brackets mean averaging over source quantum numbers 7; the overline means 
averaging over emission times and is detailed in the r.h.s. Inserting (|37D into ( p8f) yields 

P 1 (p) = f d 4 x d 4 x' e ip < x - x ' ] ^2 p^' i/>^{x) ^*,{x') , (39) 

J 7,7' 

where p = (u(p), p). Let us introduce new time and position space variables 

Y = -(x + x f ), y = x-x'. (40) 

We now define the single particle Wigner density S(X, K) of the source as 

S(Y, K) = J d'y e iK y J2 Prf ^ ( Y + I) %' ( Y ~ 1) • ( 41 ) 



7,7 



Using the hermiticity of the source density matrix p 7 y it is easily shown that S(Y, K) is real. 
Thus, we come to the expression for the single-particle spectrum, which employs the source 
function: 

P x (p) = f d 4 xS{x,p) (42) 



where the integration goes over emission times. 
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2. Two-particle quantum statistical correlations without 
final state interactions 

Let us consider a two-particle state -?/> 7 emitted by the source. Its propagation to the detector 
is governed by the Schrodinger equation 

. 7 a, b, = jj-( Xa)Xb )^ 7 (x a ,Xfc,t), (43) 

where ff(x a ,x b ) = /i(x a ) + /i(x b ). The index 7 denotes a complete set of 2-particle quantum 
numbers. (In a basis of products of two wave packets these could contain the centers X a , X^ 
of the wave packets of the two particles at their freeze-out times t a , t b , respectively.) Equation 
(f43D is solved by 

7/> 7 (x a , Xft, t) = e -^^a^)(t-to) ^ (Xq5 ^^ t()) (44) 

in terms of the two-particle wave function at some initial time t - We will assume that the 
detector measures asymptotic momentum eigenstates, i.e. that it acts by projecting the emitted 
2-particle state onto 

0™* P6 (x a , X6, t) = e*^*-*--"-*) e *(P»-*»-«6*) , (45) 

where u> 0) 6 = p„ b /2m 2 . We will only consider the case of pairs of identical particles, m a = m\, = 
m. The measured two-particle momentum amplitude is then 

^ 7 (Pa, Pft) = Hm / d 3 ^a d 3 a; 6 0°^ 5 (x a , x 6 , t) ^ 7 (x a , x 6 , t) . (46) 

If we consider identical particles, then the two-particle wave function ?/> 7 (x a ,X6, t) should be 
symmetrized (antisymmetrized). Writing this explicitly, we obtain 

A 7 (p a , Pb) = Km -= / d 3 x a d 3 x b 0°^£ 6 (x a , x&, t) [^ 7 (x a , x 6 , t) + e ia V> 7 (x fe , x a , tj\ , (47) 

where a = (a = it) for identical bosons (fermions). Relabelling the variables of integration 
in the second term on the r.h.s. of this equation as x a — > X5 and x& —>■ x a , we come to the 
expression: 

^ 7 (Pa,Pb) = lim 4= f d 3 X a d 3 X b [ e «(P--^-^0 ^(Pb^-c^t) + e ia e i(p«-xj- W „t) e i(p».x - W6 t)j* 

x V 7 (x a ,x fe ,t). (48) 

This simple algebra results in the evident conclusion that it does not matter which wave function 
should be symmetrized (antisymmetrized) at asymptotic time t: two-particle wave function 
^ 7 (x n , X6, t) or out-state wave function. We show below that the same holds for emission times. 

Using the evolution operator (see (0)), Eq. ( pf ) can be expressed in terms of the two- 
particle wave function i\) 1 emitted at earlier times as 

A 7 (p a ,p 6 ) = lim [d 3 x a d 3 x b [ e -i%»«)(*o-*)^t ( Xfl;X6) J* ^ 7 (x a ,x b ,£ ) • (49) 

t — 5-00 J ]_ 
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This is correct for all times to > max [t a , t b ], where t a j, are the freeze-out times for the two par- 
ticles. Note that the time evolution operator has been shifted (by Hermitian inversion of the 
unitary evolution operator) from the emitted two-particle state with arbitrary quantum num- 
bers 7 to the two-particle momentum eigenstate 0°"* (x a ,x&, £) which is thereby transformed 
into a plane wave at time to- It should be pointed, that such an evolution of the out-state from 
t = oo to t = to if one includes the two-particle interaction into the Hamiltonian H(x a , x&), 
would bring a distorted wave at time t instead of the plane one [|J. 

We assume that the two particles are emitted independently, implying that, at some freeze- 
out time t a , the two-particle wave function -?/> 7 (x a ,x b , t) factorizes 

^ 7 (x a ,x fe ,t a ) = — = [V 7a (x a ,t a )^ 76 (x 6 ,t a ) + e"> 7a (x fe ,t a )V 76 (x a ,t a )] , (50) 

The indices j a , , y b on the 1-particle wave functions now label complete sets of 1-particle quan- 
tum numbers. Time t a is the emission time of the latest emitted particle. At this time the 
first emitted particle of the pair has already propagated for a time t a — t b if it was emitted at 
t b < t a . During this time the first emitted particle cannot "see" the second particle as a separate 
entity, but only as part of the remaining fireball. That is why the two-particle symmetrization 
(antisymmetrization) can be effectively extracted from many-particle symmetrization (anti- 
symmetrization) after emission of the second particle. In other words a factorization of the 
two-particle wave function from the many-particle wave function, which corresponds to the 
total system, is possible only when the second emitted particle is frozen out from the fireball. 
This means that we adopt concept of the two-particle amplitude factorization. All these con- 
siderations complicate when we start with the two-particle Hamiltonian that contains terms 
where the coordinates are entagled, for instance the two-particle potential energy K(x a — x b ); 
but a discussion of this point is beyond the scope of the present letter (see j|] for details). 

Owing to the particle-particle emission symmetry and to the commutation of the free evo- 
lution operators, we can write the measured two-particle momentum amplitude as 



A ia,i b (P<^ Pb, t a , h) = lim — = / d 3 x a d 3 x b (p°^* b (x , x 6 , t) 



x r-i[h(x a )(t-t a ) + h(x b )(t-t b )] , , t s , , t s + 



7a V^ffl' v aJ Yj b * 

7a V X &' ta) Y~y b \ 



+ e ia e -i[h(x b )(t-t a )+h( Xa )(t-t b )] - (x ^ ta) . (Xa; tb)] _ (51) 



Now, we first relabel the variables of integration in the second term on the r.h.s. of this equation, 
what results in symmetrization (antisymmetrization) of the out-state. Secondly, we invert the 
evolution operator, which therefore acts on the out-state and brings it from asymptotic time 
t to initial times t a and tb- As a result, we come to the final expression of the two-particle 
amplitude: 

A la , lb (p a , Pb ,X a ,X b ) = -j= J d 3 X a d 3 X b [ e i(Pa-Xa^ b -X h ) +e -4a e i(pa^^ b -Xa)j ^ a ( Xa ) ^{x b ) , 

(52) 

where x° a = t a and x° = t b . We therefore represent the measured two-particle momentum 
amplitude as a projection of a non-symmetrized two-particle wave function taken at emission 
times onto symmetrized (antisymmetrized) plane waves taken at emission times as well. 

16 



The two-particle probability is obtained by averaging (|5^) and its complex conjugate with 
the density matrix defining the source. This density matrix is characterized by a probability 
distribution for the two-particle quantum numbers (j a , j b ) and by a distribution of emission 
times (t a , t b ). We write: 



^2(Pa,Pfe) = (|A 7ai7b (p a ,p 6 ;£ a ,£ 6 )| 2 ) = (53) 

\ ' "la ,1b 

= Z_^ I dta dU dt a i dty p 7a7a , p 7b7b , ^4 7a , 7b , (Pa, P&; t a ',t b i) A* lalb (p a , p b , t a , t b ) , 

7a7b,7 a '7(/ 

where anglular brackets on the r.h.s. of this equation mean averaging over the quantum numbers 
7 a , 7;, and the overline means averaging over the initial (emitting) times. We made the ansatz 
Pia-y b n a ny = Piai a i Pibi b n which factorizes the initial density matrix p 7a7b , 7a , 7b , in such a way that 
independent emission of the two particles is ensured. 

According to fl52f) , the probability consists of four terms, which we write as 

P 2 (p a , Pft ) = P n + P 22 + P 12 + P 21 , (54) 

which have the structure (a\ + a 2 )(ai + a 2 )* = aia* + a 2 a-2 + aia^ + a 2 a:j\ First diagonal term 
reads 



x j d A x b dS b e i{pb - Xb ^ yb) (^ lb {x b )r^)) ■ (55) 



We introduce the new momentum variables if and p as 

1 p p 

K =^(Pa+Pb), P = Pa~Pb => p a = K+-, p b = K - - . (56) 



Using these variables and the space-time variables (fTOfl, P n can be rewritten in the form 



2/ / -..„ 



Pn(p, K) = \Jd A X J d*x e< K +y>~ (tf> 7o (x + |) < (* - 

x |rf 4 r /Ae^-^(^ 76 (r + |)^;(r-|))^ . (57) 

In the integrals over the variables a; and y, we immediately recognize the freeze-out Wigner 
density (0). But before rewriting this term let us note that P 22 can be obtained from Pn by 
a mutual change of momenta p a ^ p b , which in turn results in changing p to —p. Hence, to 
obtain P 22 from Pn we need only to change the sign before p. Performing such a change, we 
come to an equality of the diagonal terms Pu(p,K) = P 22 (p, K). Putting together these two 



terms, we obtain 



P n (p,K)+P 22 (p,K) = j ' d i Xs(x,K + 1P ^ jd 4 YS(Y,K-^j , (5* 
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i.e. the product of two single-particle probabilities (see Eq. (|42"|)) to registrate independently 
two particles with asymptotic momenta p a = K + p/2 and p^ = K — p/2, respectively. 

We now turn to considering of the cross term P 12 = aia^j, which is a complex conjugate to 
the second cross-contribution P 2 i = a.iQ-2, hence their sum is real. We have 



x f d 4 x b d*y b e^ x >- p ^ (iJ lb (x b )rjy b )) ■ (59) 



By the same change of variables as for diagonal terms, the exponents in the integrand can be 
rewritten as: (p a ■ x a - p b ■ y a ) = (p ■ X + K ■ x) and (p b -x b -p a - y b ) = (-p -Y + K -y). Using 
the definition of the source function (Ell), we obtain: 



p 12 (p, R) = ]- e ia ! d A X e ipX S(X, K) f d A Y e~ ip - Y S(Y, K) . (60) 

The second cross term P21, which is complex conjugate to P12, is proportional to the exponent 
exp (— ia), but the remaining factor is real. Hence, the sum of two cross terms looks like 
P12 + P21 = \Ae ta + \Ae~ %a = A cos a, where A is real. We are now ready to write the total 
expression for the two-particle probability 

P 2 (p,K)= Jd 4 Xs(x,K + ^j Jd A YS(Y,K-^j 

+ cosa f d 4 Xe ip - x S(X,K) f d 4 Ye- ip - Y S(Y,K) , (61) 

and we finally get the two-particle correlator as 

\[d A Xe ip - x S(X,K)\ 2 , x 

C(p, K) = 1 + cos a ^ r v ' ' r- , (62) 

jd^XS[X,K + l) jd 4 YS[Y,K-lj 

where the source function S(X, K) is defined in accordance with Eq. (|4l]) and all integrations 
are taken at emission times or on the freeze-out hyper-surface. 
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